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There exists a class of gauge models incorporating a finite density of matter in which the Higgs 
mechanism is provided by condensates of gauge (or gauge plus scalar) fields, i.e., there are vector 
condensates in this case. We describe vortex solutions in the simplest model in this class, the 
gauged 517(2) x U(1)y u-model with the chemical potential for hypercharge Y, in which the gauge 
symmetry is completely broken. It is shown that there are three types of topologically stable vortices 
in the model, connected either with photon field or hypercharge gauge field, or with both of them. 
Explicit vortex solutions are numerically found and their energy per unit length are calculated. The 
relevance of these solutions for the gluonic phase in the dense two-flavor QCD is discussed. 

PACS numbers: 11.27.+d, 11.15.Ex, 11.30.Qc, 12.38.-t 


I. INTRODUCTION 

Vortices are very special soliton-like excitations. In condensed matter physics, vortex solutions were discovered 
by Abrikosov pj: They play a crucial role in the dynamics of type II superconductors. In relativistic quantum field 
theory, vortices were first considered by Nielsen and Olesen ||. They can be important in cosmology, astroparticle 
physics, and for the confinement dynamics in QCD-like theories (for reviews, see Ref. 0). Vortex and vortex-like 
solutions in relativistic models at finite density of quark matter were studied in Ref. dj. 

The existence of vortex solutions is usually connected with the Higgs mechanism in abelian (7(1) gauge dynamics. 
The major degrees of freedom are a complex scalar Higgs field and a (7(1) gauge field. In this paper, we will describe 
vortex solutions in a somewhat different class of gauge models incorporating a finite density of matter. The Higgs 
mechanism in these models is provided by condensates of gauge (or gauge plus scalar) fields, i.e., there are vector 
condensates in this case. This class in particular includes the SU(2)l x U(1)y electroweak theory in the presence 
of a superdense fermionic matter 0, @j, (nonrenormalizable) models including massive vector bosons with a chemical 
potential for electric charge 7|, and the gauged linear SU(2)l x (7(1 )y u-model (without fermions) with a chemical 
potential for hypercharge Y [|. 1 

It is also noticeable that a recently revealed gluonic phase in neutral two-flavor quark matter nSf also relates to 
this class. In the gluonic phase, vector condensates of gluons cure a chromomagnetic instability |l3| in the two-flavor 
superconducting (2SC) solution and lead to spontaneous breakdown of the SU(2) C x U(l) em x SO( 3) ro t symmetry 
down to SO( 2) rot . Here SU(2) C and U(l) em are the color and electromagnetic gauge symmetries in the 2SC medium, 
and <SO(3) ro t is the rotational group (recall that in the 2SC solution the color 5(7(3) c is broken down to SU(2) C ). In 
other words, the gluonic phase describes an anisotropic medium in which the color and electric superconductivities 
coexist. Because it is naturally to expect that cold quark matter may exist in the core of compact stars (for reviews, 
see Ref. |h|), it would be of great interest to describe the dynamics of this phase in more detail, in particular, to 
clarify whether vortex excitations exist there. 

Since the dynamics of the gluonic phase is very rich and complicated, as a first step, it would be appropriate to use a 
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simpler model relating to the same universality class, i.e., with the same sample of the symmetry breaking. Fortunately, 
such a model exists: It is the gauged a-model with the chemical potential for the hypercharge introduced and studied in 
Ref. @j. Let us describe it in more detail. Its Lagrangian density reads (we use the metric g = diag(l, —1, —1, — 1)): 

C = -^F^F^ a) - ^F^F^ y) + [(£>„ - ifiSuo)$]HD v - */A°)$ - - A($ t $) 2 , (1) 

where the covariant derivative — igA^ — (ig'/2)B IJ ,, $ is a complex doublet field <I> T = (<p + ,tpo)i and the 

chemical potential g is provided by external conditions (to be specific, we take g > 0). Here A p = A“r°/2 are 5/7(2)^ 

gauge fields (r a are three Pauli matrices) and the field strength F$ = d^A^ — d v A^ + ge abc A^A^. B M is a Uy(l) 
gauge field with the field strength F^ = d^B^ — d v B^. The hypercharge of the doublet <I> equals +1. This model has 
the same structure as the electroweak theory without fermions and with chemical potential for hypercharge Y. Note 
that the terms with the chemical potential are SU(2)l x U{1)y (and not SU(2)l x SU(2)r) symmetric. This follows 
from the fact that the hypercharge generator Y is Y = 2/|, where 1^ is the third component of the right handed 
isospin generator. Henceforth we will omit the subscripts L and R, allowing various interpretations of the 5/7(2). 

Because the chemical potential explicitly breaks the Lorentz symmetry, the symmetry of the model is SU (2) x 
U(1)y x SO(3) rot . As was shown in Ref. J§j, for sufficiently large values of the chemical potential g , the condensates 
of both the scalar doublet $ and the gauge field A^ occur. These condensates break the symmetry 5/7(2) x J7(l)v x 
50(3) ro t down to 50(2) rot . This is the same symmetry breaking pattern as in the gluonic phase of dense QCD jl^ . 
while the gauged tx-model is much simpler: It is renormalizable and for small coupling constants g , g' and A, the tree 
approximation is reliable there. 

As will be shown below, topologically stable vortex solutions exist in this model indeed. Their structure is much 
richer and more complicated than that in the Abelian gauge model. In particular, there are different types of vortices 
connected either with photon field or hypercharge gauge field B or with both of them. 

The paper is organized as follows. In SectionEl a general analysis of vortex solutions in this model is realized. We 
discuss a topological stability of these solutions and introduce three types of vortices: Magnetic vortices, hypermagnetic 
vortices, and hybrid ones. They are connected with photon field, hypercharge field B and both of them, respectively. 
In Section EH an explicit solution for magnetic vortices is obtained. In Sections m and El explicit solutions for 
hypermagnetic vortices and hybrid ones, respectively, are derived. In Section eh we summarize the main results of 
the paper. 


II. VORTEX SOLUTIONS: GENERAL CONSIDERATION 

The equations of motion in the gauged cr-model at finite density follow from Lagrangian 0: 

-(A, - ig5vo){D u - ig6 v0 )<f> - m 2 $ - 2A(<L t $)$ = 0 , 


flV 


&*F$ + ig ( d v - i 9 -B v i^B v ) $) —t> 


, 5 _ 
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g£ abc A^b) F (o) + t$ + 2gg8 v0 & T —§ = 0 , 


( 2 ) 


( 3 ) 


d»F^ + ^d>t(A - ig6 „ 0 )$ - ^[(A - ig5 v0 m^ = 0. (4) 

As was shown in Ref. §|, for sufficiently large values of the chemical potential g, the vacuum solution of 0-0 is 
given by 




4 (3) _ 

" 71 ’ 


Wj-> = (Wj +) )* = c 


$ T = (0, v 0 ), 


( 5 ) 
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where 


sj ( g 2 + 64A)/u 2 — 8(8A — g 2 )m 2 — 3 gp 
V2(8A - g 2 ) 


( 6 ) 


Wf> = i *7^), < h T = (<A*7 <Po) and all other fields are equal to zero. 2 It is clear that this solution implies 

that the initial symmetry 5/7(2) x U( l)y x 50(3) ro t is spontaneously broken down to 50(2) rot . In particular, the 
electromagnetic U(l) em , with electric charge Q e m = I 3 +Y/ 2 , is spontaneously broken by the condensate of W bosons, 
i.e., electric superconductivity takes place in this medium. 

Note that because the U(l)y symmetry is local, for nonzero chemical potential p one should introduce a source 
term BqJq in Lagrangian density CD in order to make the system neutral with respect to hypercharge Y. This 
is necessary since otherwise in such a system thermodynamic equilibrium could not be established. The value of 
the background hypercharge density Jo (representing very heavy particles) is determined from the requirement that 
B 0 = 0 is a solution of the equation of motion for B 0 (Gauss’s law) 0, fl5| . There exists an alternative description 
of this dynamics in which a background hypercharge density Jo is considered as a free parameter and p is taken to 
be zero. Then Gauss’s law will define the vacuum expectation value (Bo). It is not difficult to check that these two 
approaches are equivalent if the chemical potential ll in the first approach is taken to be equal to the value ^(Bq) 
from the second one. In this paper, following Refs. am , we use the first approach. 

As was shown in Ref. 0 ], in accordance with the sample of spontaneous breakdown of the global symmetry, 
SO(3)rot —► SO{2) rot , there are two gapless NG bosons in this model. The other excitations are massive (the Higgs 
mechanism). Because electric superconductivity is realized in the model, it is naturally to expect that vortices may 
exist there. As will be shown in this paper, this is indeed the case. 

Vortices are topologically nontrivial configurations with fields approaching their vacuum values at spatial infinity 
)3]. Mathematically, vortices are topologically stable whenever the first homotopy group T<\(G lnt /H int ) of the vacuum 
manifold G, nt / H int is nontrivial (here Gi nt is an internal symmetry group of the action in a model and its subgroup 
Hi n t is a symmetry group of the vacuum.) We emphasize that this criterion of the topological stability is sufficient 
but not necessary. In particular, as will be shown below, magnetic vortices, the most interesting vortex solutions in 
the present model, are topologically stable due to another criterion. 

According to Eq. ©, two complex fields wj "* and po have nonzero vacuum expectation values in the ground state 
of the gauged er-modcl at finite density. Therefore a priori there can be different types of topologically stable vortex 
solutions: Both the field wj and the field po can wind around the spatial circle at infinity. I.e., in the cylindrical 
coordinates (t, as p —> oo, the asymptotics of these fields take the form 


wj-) 


Ce^, 


Vo 


v 0 e 


incf) 


(7) 


where l and n are integer. Let us turn to the analysis of a topological stability of solutions with these asymptotics. 

In this model, the internal group Gi n t is SU( 2) x /7(l)y and the symmetry group Hi nt of vacuum © is trivial. 
Therefore the first homotopy group of the vacuum manifold in the model is 7Ti (Gi n t/Hi nt ) = 7Ti (5/7(2) x Z7(l)y) = 
7Ti(5C/(2)) + 7Ti(t/(l)y) = 7Ti(U(1)y) = Z, where Z is the group of integer numbers (here we use the fact that since 
the 5/7(2) manifold is simply connected, 7Ti(51/(2)) = 0). This seems to suggest that there exist only topologically 
stable vortex solutions connected with the hypercharge Z7(l)y. However, as will be shown below, the situation in 
this model is more sophisticated and there exist also other topologically stable vortices. The point is that while the 
scalar field d> is assigned to the fundamental representation of the 517(2) group, the gauge field A^ is assigned to 
the adjoint one. Because of that, 5/7(2) gauge transformations for A^ reduce to 50(3) ones. It is important that, 
unlike 5/7(2), the manifold of 50(3) is not simply connected and 7Ti(50(3)) = Z 2 with Zi = 0,1. As will become 
clear in a moment, this point is quite important in the analysis of the topological stability of the solutions. 


2 Here “sufficiently large values of means the following: When m 2 > 0. //. should be larger than the critical value // rr = m, and for 
m 2 < 0, // should be larger than ji, rr = fj\rn\/2\fX (the critical value g\rn\/2\/\ coincides with the mass of W boson in the vacuum 
theory with fi = 0 and m 2 < 0). 
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We begin by considering the conventional case when l = 0 in asymptotics 0 : 

- C, v?o - v 0 e in *. (8) 

This configuration can be obtained from the vacuum configuration © by using the U{\)y gauge transformation e m ^ . 
It is clear that in this case the solutions with different n are assigned to different topological classes described by the 
homotopy group 7Ti(t/(l)y) = Z. Therefore, if such vortex solutions exist, they are topologically stable. Due to the 
evident reasons, we will call them hypermagnetic vortices. They possess hypermagnetic fluxes with a (hypermagnetic) 
winding number H = n. 

Let us now turn to a more interesting case of asymptotics 0 with n = 0: 

-> Ce"*, ipo^vo. (9) 


It is a rather unusual situation: The phase lip is now related to a vector field and not to a scalar one. It is clear that 
configuration © can be obtained from vacuum configuration © by using the gauge transformation 

e -ii0Q em _ e -i/0T 3 /2 e ~il(t>Y/2 


where electric charge Q em = I 3 + Y/2 is equal to -1 for the field and zero for the scalar field ipo- 

Because for the vector fields the SU (2) transformations are reduced to the 50(3) ones, it is important to distinguish 
two cases: with l being even and l being odd. In the first case, with l = 2k, as the angle <f> running from 0 to 
27t, the gauge transformation m defines closed loops both on the manifold of 51/(2) and that of U(l)y. Then, 
because 7Ti (517(2) x U(1)y ) = ki(U(1)y) =Z, we conclude that these solutions are topologically equivalent to the 
hypermagnetic vortex solutions with the winding number H = —k considered above. 

The case with odd l = 2fc+l is more sophisticated. In this case, as (j> changing from 0 to 27r, the gauge transformation 
03 does not define closed loops neither on the SU(2) manifold nor on the U(1)y one. Therefore in this case one 
cannot use the homotopic group 7Ti (51/(2) x U(I)y) for studying the topological stability of these solutions. Let 
us show that for this purpose one can use the homotopic group 7Ti(50(3)) instead. The reasons are the following. 
Because the hypercharge of the gauge field A^ is zero, the action of the gauge transformation 03 on it is given by 
the operator e _ ^ 2fc+1 ^ / 3 J / 2 , where I'p dj is a SU( 2 ) generator in the adjoint representation. As the angle (p changing 
from 0 to 2-7T, this transformation does define a closed loop on the 50(3) manifold, which coincides with the vacuum 
manifold of . 3 Therefore, one can indeed use the homotopic group 7Ti(50(3)) = Z 2 for the description of the 
topological stability of these solutions. 

Now, because odd values of l correspond to the element +1 of 7Ti(50(3)) = Z 2 , we conclude that this loop cannot be 
unwinded and, therefore, these vortex solutions are topologically stable. Due to the form of gauge transformation 03 . 
we will call them magnetic vortices. They possess magnetic fluxes with a (magnetic) winding number M = — (2k + 1 ): 
The minus sign here corresponds to a negative electric charge of W~ bosons. A detailed study of the topological 
structure of the set of these solutions is beyond the scope of this paper. Here we will only argue that the solutions 
with different values of k are topologically inequivalent. The argument is based on gauge transformation 111 Oil . It 
implies that although for odd l = 2 k+1 the path defined by this transformation on the U(l)y manifold is not closed, 
it contains k windings. Therefore the solutions with different k should be assigned to different topological classes, i.e., 
they cannot be transformed from one into another by using a continuous deformation. 

The last case is that with the “hybrid” asymptotics 

Wj~ } -► Ceipo -> v 0 e in<1, (11) 

with both l and n being nonzero. In this case the gauge transformation 113 is replaced by 

e -il<pQem & in<pY _ e ~il<t>T 3 /2 e i(2n-l)</>Y/2 Q 2 ) 


3 These two manifolds coincide because vacuum configuration © includes two linearly independent isovectors, A ^, and , and there 
is no subgroup in SO (3) under which two linearly independent isovectors are invariant. 
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Using the same arguments as above, we conclude the following. For even l = 2 k, these solutions are topologically 
equivalent to the hypermagnetic vortices with fluxes corresponding to the winding number H = n — k. As to odd 
l = 2k + 1, they yield new topological classes of solutions. These solutions carry both magnetic and hypermagnetic 
fluxes with winding numbers ( M,H ), with M taking odd integers —l = — (2k + 1) and H = n. It will be shown in 
Sec. |V| that the hybrid vortices can be considered as composites of the magnetic and hypermagnetic ones. 

As to the topological structure of the solutions for hybrid vortices with odd l = 2k + 1, a comparison of Eq. m 
with Eq. ITUli suggests the following generalization of the conclusion that was made above for the case of the magnetic 
vortices: The hybrid vortices with different n — k are topologically inequivalent. 

This concludes the analysis of the topological stability of vortex solutions. As will be shown below, these three 
types of vortex solutions exist in the present model indeed. 


III. MAGNETIC VORTICES 


In this section, we will analyze the solutions for the magnetic vortices having asymptotics © with l being odd. To 
be concrete, we will take l = — 1. In accordance with the discussion in the previous section, one should expect that 
the magnetic winding number M equals +1 in this case. 

Vortex solutions are static ^-independent configurations. We will consider the following ansatz for them in the 
cylindrical coordinates (i, p, <fi, z) (the reasons for choosing this ansatz will be discussed in more detail below): 


Wi \t,p,</>,z) = 


W(p)e~ 

~7T 


A t 3) (t, P, <t>, z) = A { t 3) (p), A^t, p, (j), z ) = A^(p), 

p, <j>, z) = Z^p), <p 0 (t, p, <j>, z) = v{p), 


(13) 


g' ^( 3 ) -\-gB 

where Ap and Z$ are angular components of the electromagnetic field = —^== ^ and the Z-boson field 

Z^ = aA ^j ^ (all other components of the fields are equal to zero). 4 

Using equations of motion ©-©, one can check that ansatz da is consistent with them and they reduce to the 
following system of five coupled nonlinear ordinary differential equations: 


W 


W_W 2gWA t +gZ t )W 2 / <„ )2 _ W^+gZ^ _ £ ^ = 

p p V 9 2 + g 2 ) 2 


(14) 


(Af)" + - 9 2 W 2 A f - V + gpv 2 = 0 , 

p 2 


(15) 


A n + _^± gg ' w2 M _ gW A <t> + gZ<t,) \ = n 

* p p 2 AF+^pyp v 9 r + 7 I ) ’ 


(16) 


ry,, , Z ’<t> Z 0 g 2 +g' 2 r7 2| g 2 w 2 (l g{g'A <j> + gZ <k )\_ n 

* — z * v 7FW r 


(17) 


// , 9 2 +g 2 ry2„, 5 2 w 2 


v H- 

p 


4 + 


g A 


( 3 ) 


— \i I v — mrv — 2Xv 3 = 0 . 


(18) 


4 Recall that in the cylindrical coordinates, a vector field A = (A x , A y , A z ) is decomposed as A = A p p + A^cf) + A z z, where A p 
A x cos 4> + A y sin <j>, A $ = —A x sin (f> + A y cos 4>. 
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Note that although the phase in the exponent e *7 in wi ^ in Eq. m is of the electromagnetic U(l) em origin, this 
ansatz contains in addition also Z^ field. This is because the W condensate mixes the gauge fields A M and Z M (see 
Eqs. Hflll and 1171) '). 

It is instructive to compare ansatz m and equations of motion m-m with those in the Abelian Higgs model. 
Obviously, W i - '* plays the role of the Higgs field. The fields A$ and Z$ are analogous to the Abelian gauge field. The 
fields A [ 3 -* and v do not have analogs in the Abelian Higgs model. Of course, they are present in ansatz m in order 

(o\ 

to ensure the correct vacuum asymptotics at infinity. Thus, the fields A\ ; and v have a somewhat different status in 
the vortex ansatz comparing to the fields Wj \ A^, and Z^. We will discuss this point in more detail below, when 
we consider the ultraviolet boundary conditions for these fields. 

The infrared boundary conditions at p —> oo follow from the requirement of finiteness of the energy per unit length 
of a vortex, which implies that fields tend to their vacuum values: 

W ( p) -► - -y , ^* 3) (p) ~^= i v{p)->v o, A 0 (p)-► 0, ^(p)->• 0. (19) 

As to the ultraviolet boundary conditions at p = 0, because p = 0 is a regular singular point of the system of 
differential equations m-m, these boundary conditions follow from the equations themselves if one utilizes the 
condition of regularity of the solution at this point: 

W(0) = A 0 (O) = Z 0 (O) = 0, (20) 


(4 3) )'(0)=</(0)=0. 


( 21 ) 


The ultraviolet boundary conditions lt20l) for the fields W and A^, Z$ are like those for a Higgs field and a gauge field, 
respectively, in the Abelian Higgs model: they guarantee that these fields are single valued at p = 0. Unlike them, 

(o\ 

the fields A) ; and v have Neumann boundary conditions which imply that there is no influx of momentum for these 
fields at the boundary p = 0. 

We solved numerically Eqs. m-m with boundary conditions tun-® by using MATLAB. To be concrete, we 
considered the case with m 2 > 0 and chose the parameters p and gv o /v/2 as in Ref. ||: p/m = 1.1 and gvo/\/2m = 0.1. 
For the values of the coupling constants g and g' we used those from the electroweak theory: g = 0.65 and g' = 0.35. 
Then, using Eq. ©, one gets the value of the coupling constant A: A = 0.53. Therefore we work in the weak coupling 
regime when the tree approximation is reliable. 

(ON 

The solution is shown in Fig. 1, where dimensionless fields W/m , A\ /m , A^/m, Z^/m, v/m are plotted against 
r = pm. The boundary conditions m were set at r max = 80 and it was checked that starting from r max = 40 the 
form of the solution practically does not depend on the choice of 7'max- 

Substituting ansatz m into Eq. Q. we calculate the effective potential (V = — C): 


(W') 2 W 2 (l gjg'A^ + gZ^ Y (Ap) 2 
2 2 (p vV + 9 2 ) 2 


g^fW 2 , 1 ( A , 

2 + 2 V 4 * 



+K) 2 + 


r< 2 2 
'9 ry 2 2 , _2 2 , \ 4 , 9 Tir2 2 
■ - ZiV +771 V + AV +— W V . 



( 22 ) 


Subtracting the energy density of the vacuum solution, we find that the energy of the magnetic vortex per unit length 
is equal to 1.27 in units of m 2 . 

Note that the 517(2) x U(1)y invariant definition of the electromagnetic field strength is 


P (em) _ 9F& 2 g' 

^ v / 7+F ^fg ^ W i $t<i> 


( 23 ) 
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FIG. 1: Dimensionless fields 


a (3) 


-, and — as functions of r = pm for the magnetic vortex solution. 


For the vortex solution with ansatz CC9. it yields the conventional relation for the magnetic field H: 

H = curl A . 


(24) 


It is well known that the magnetic flux of vortices is quantized and related to the winding number, which is a 
topological invariant. It is easy to check that this property takes place in the present case. Indeed, as follows from 


Eqs. itTTfl) and G3b while the field A^ has the asymptotics A- as p —> oo (with e = 


aa 


=), the field Z$ rapidly 


e P ' ' \/g 2 +g'- 

decreases in this limit. Therefore we find that for the magnetic vortex only the flux related to the electromagnetic 
field A^ is nonzero. Then, using Eq. in , we get: 


= J dxdyH z = j> dx ■ A = 

I.e., as was expected, this vortex solution corresponds to the winding number +1. 


. i 27r 

cLcp P-A(p | p=oo — • 


(25) 


IV. HYPERMAGNETIC VORTICES 


A hypermagnetic vortex possesses asymptotics ©• To be concrete, we will take n = — 1. As follows from the 
discussion in Sec. m the hypermagnetic winding number H should be now equal to -1. 

In this case, we will use the following ansatz: 

Wt ) (t,p,<t),z) = -^ 7 =-^, (t, p, <f>, z) = a[ 3) (p), = A^p), 

Z<t>(t, p, (j>, z) = Z^p), y 0 (t,p,<t>,z) = v{p)e~ l<1, . (26) 

For this ansatz, the equations of motion ©-Cl yield 

W" + w + g 2 ( {A f)f {g'At + g z+y \ w _ £ Wv2 = 0 

p V g 2 + g 2 J 2 


( 27 ) 















(28) 


t(3)y/ , (^t 3) )' ^2 tt/2 a (3) 9 2 A { 3)„,2 


«T + 


-^T-t^v + s ^ = o 1 
P 2 


^4" _|_ (ff ~^0 + _ Q 


01 p P 2 


g 2 +9 


(29) 


■^0 ^0 . V 9 2 +; 


r^// cp — <-p 

^ 4> ' 2 

p p 


- 0 - n/T+T^) 


t; 2 - 


g 3 (g'^0 + g-^>) 

g 2 + g' 2 


it 2 = o, 


(30) 


7/ 1 / 2 

P 4 


Q - \/g 2 + g' 2 0) 




— /u I v — m 2 v — 2Xv 3 = 0. 


(31) 


The infrared boundary conditions for the hypermagnetic vortex ansatz at p —> oo are the same as for the magnetic 
vortex and given in Eq. (HU . The ultraviolet boundary conditions at p = 0 are 


,(0) = A^O) = Z^(0) = 0, 


(32) 


(4 3 ) )'(0) = fr(0) = 0. 


(33) 


Comparing these boundary conditions with those in Eqs. m and (ZD for the magnetic vortex, one can see that 
the ultraviolet boundary conditions for W and v are now interchanged. Of course, this point is related to the fact 
that the phase in the scalar field tpo (and not vector field W) is relevant for the hypermagnetic U( l)y vortex. The 
numerical solution of Eqs. |77 |) -P 1) with boundary conditions m cub and iCTl is shown in Fig. 2, where we use 
the same set of parameters as for the magnetic vortex. 

The effective potential for ansatz (El is 


V = 


(W') 2 1 g 2 (g'Ar/) + gZ<t ,) 2 _ w 2 _ (Af) 2 


g- +g- 


g 2 (A ( V) 2 W 2 1 / , A 


2 l^ + T 


-Jz:, + ^\ - 


gA 


(3) 


P v 


+(v') 2 + i 0 - \Jg 2 + g' 2 Z^j v 2 + m 2 v 2 + Xv 4 + ^W 2 v 2 . 


(34) 


Subtracting the energy density of the vacuum solution, we find that the energy density of the hypermagnetic vortex 
per unit length in units of m 2 is equal to 0.79. 

It follows from Eqs. El and m that A ^ and Z^ have the asymptotics- 2g and — . 2 as p —> oo. 

g'p\/g 2 +g 12 p\/g 2 +g' 2 


Therefore, the fluxes of both fields A^, and Z$ are nonzero 


f(H) 

u pm. 


47i -g 

g'Vg 2 + g ' 2 ’ 


(35) 


t (h) _ 

T Z — 


47r 

vV + g ' 2 ' 


( 36 ) 
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FIG. 2: Dimensionless fields —, 

m ’ 



m 



m 5 



m ’ 


and 


as functions of r = pm for the hypermagnetic vortex solution. 


It is not difficult to see that these two fluxes correspond to the hypermagnetic flux of the U(l)y gauge field (the 

/'o't 

hypermagnetic flux of the A^ ; field is zero): 


{H ) _ _4tt 

r - y 


(37) 


Note that because the coupling of the Higgs field with the field B^ in the covariant derivative D M — igA^ — 
(ig 1 /2)B I1 is expressed through g '/2 rather than g ', this flux corresponds to the winding number -1, as was expected. 


V. HYBRID VORTICES 

The hybrid vortex solutions carry both the magnetic and hypermagnetic fluxes with winding numbers (M, H). Their 
asymptotics are given in Eq. CD- To be concrete, we will consider the solutions with l = —1 ,n= qH corresponding 
to (M, H) = (1, =pl), i.e., we consider hybrid vortices with the equal and opposite phases of the Wz~ 1 and tpo fields. 
Clearly, the ansatz for these vortices should have the form 

Wz~\t, p, <j>, z) = W<yP ^ -> A 3) (t,P><l>,z) = ^ 3) (P), A^pifaz) = A^p), 

Zyp^,z) = Z^p), <p 0 (t,p,<j>,z) = v^e-^ , (38) 


where <r = ±. 

The equations of motion for these vortices are sort of a combination of the equations of motion for the magnetic 
and hypermagnetic ones: 


W h 


W' 

p 


w 

T2 


Vg 2 + 9 2 


+ 9 ( A t ) - 


(3)n2 + gZ^Y 


9+3 


W - —Wv 2 = 0, 
2 


P 


(39) 
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(4 3) )" + — - s ! »- 2 4 3) - +»» 2 = 0 ■ 

P 2 


(40) 


A'L ^ ^ - - 


+ gg'^ 2 f i _ g(g'7 + g7) ] = 0 


p p 2 vV + g' 2 \ p vV + g' 2 J 


(41) 




^ - % + yg+g ( n _ +-^7 fi - TV 7) 

p p 2 p vV+7 2 \ p vV + g 2 / 


= 0 , 


(42) 


V P 4 


v r , 2 (a \ 

- l[W 2v + ( —^- Pj v ~ m2y — 2Au 3 = 0 . 


(43) 


The infrared boundary conditions at infinity remain the same as before. Requiring the regularity of solutions at 
p = 0 leads to the following ultraviolet boundary conditions: 


W(0)=«(0)=A*(0) = Z*(0) = 0, 


(44) 


( 4 3 ) )'( 0 ) = 0 . 


(45) 


Note that in contrast to the cases of the magnetic and hypermagnetic vortices, now both functions W(p) and v(p) have 
the Dirichlet boundary conditions at p = 0. The reason of this is clear: Because both w7 and <po have nontrivial 
phases for hybrid vortices the regularity of solutions at p = 0 can be ensured only if the moduli of these fields 
vanish at this point. 

The numerical solutions of Eqs. with boundary conditions m , < eu >, and (TT51i are shown in Fig. 3 and 

Fig. 4. For these solutions, the same set of parameters was used as in the previous cases. 

The effective potential for ansatz EH is equal to 


(W'f W* 1 _ gjg'A^ + gZt) 2 _ (A {3y ) 2 

2 2 v / 7T7^ 2 


g 2 {A ( t 3) ) 2 W 2 

2 




+(v’) 2 


1 / 2er 

4 V7 


2 2 
V 9 2 + g' 2 Z<t> ) ^ 2 + iti 2 v 2 + Av 4 + ^W 2 v 2 . 


(46) 


Subtracting the energy density of the vacuum solution, we find that in units of m , the energy per unit length of 
hybrid vortices with the equal and opposite phases are equal to 1.49 and 2.3, respectively. 

It follows from Eqs. (1411) and 1421) that Aj, and Za, have the asymptotics 9 ( 1 ~ 2cr ')+ 9 and — . 2a as p —> oo. 

pgg'\/9 2 +9' 2 psjg 2 +g' 2 


Therefore, the fluxes of hybrid vortices are equal to 


■p(MH) 

‘'em 


2tt [g 2 (1 — 2cr) + g' 2 ] 

gg'Vg 2 + g' 2 


(47) 


T (MH) 

r z 


Ana 

Vg 2 + g' 2 ' 


( 48 ) 
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FIG. 3: Dimensionless fields —, 

m 1 



m 



m ’ 



m 5 


and 


V 

m 


as functions of r = pm for the hybrid vortex solution with equal phases. 



(3) ^ ^ 

FIG. 4: Dimensionless fields —. ——, —, —, and — as functions of r = pm. for the hybrid vortex solution with opposite 

m 1 m ’ m ’ m 1 m ' J 1 1 

phases. 
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These expressions can be rewritten in the following transparent way: 

'p(MH) _ i tt(M) i 

u em em 1 em j 

(49) 

t { z 1h) = 1 • f { z m) + a ■ b { z h) , 

(50) 


where and z are the fluxes of the magnetic and hypermagnetic vortices given by Eqs. iHQ) and i35i) . inst - 

respectively (recall that = 0). This suggests that the hybrid vortices can be considered as composites of the 

magnetic and hypermagnetic ones. Note that the sum of the energies per unit length for the magnetic vortex and the 
hypermagnetic one is equal to 1.27 + 0.79 = 2.06. On the other hand, because the energies for the hybrid vortices 
with a = + and er = — are equal to 1.49 and 2.3, respectively, they satisfy the inequality 1.49 < 2.06 < 2.3. This 
inequality implies that while the hybrid vortex with equal phases is stable, the hybrid vortex with opposite phases is 
not. 


VI. CONCLUSION 

The set of vortex solutions we obtained in the gauged SU( 2) x U(1)y er-model with the chemical potential for 
hypercharge is quite rich. In particular, there are different types of vortices connected either with photon field 
or hypercharge gauge field, or with both of them. The richness of this set is provided by the structure of the 
ground state which includes vector condensates. This ground state describes an anisotropic medium with electric 
superconductivity. 

It is quite noticeable that the sample of symmetry breaking in this model is the same as in the gluonic phase in 
neutral two-flavor QCD PHI . while the present model is much simpler. The point is that because of a nonzero baryon 
density in the gluonic phase, its dynamics in the hard-dense-loop approximation is mostly provided by quark loops. 
And since most of the initial symmetries (including rotational SO(3) ro t ) in this phase are spontaneously broken, the 
calculation of its effective action is quite involved. This makes the study of vortex solutions in the gluonic phase to 
be quite a difficult problem. On the other hand, vortices could be very important for understanding the dynamics of 
quark matter in compact stars. The existence of such solutions in a related but much simpler model is encouraging. 

It is instructive to compare the present vortex solutions with those in such complex condensed matter systems as 
the B phase in 3 He [01 and high T c superconductors E2- Unlike the conventional Abrikosov-Nielsen-Olesen vortices 
Qi, the initial symmetries are not completely restored in the core of vortices in those systems. In 3 He, vortices 
with an A-phase core or with asymmetric core are realized. In high T c superconductors, there is an antiferromagnetic 
condensate in the vortex core. A similar situation takes place for the magnetic and hypermagnetic vortices in the 
present model. In the core of a magnetic vortex, while the order parameter wi ^ is zero, the order parameter ipo is 
not (see Fig. 1). Therefore, in its core only the U(l) em symmetry is restored. In the core of a hypermagnetic vortex, 
we found that while (fo = 0, the order parameter Wz 1 is not zero (see Fig. 2). Therefore, only the U(1)y symmetry 
is restored in the core in this case. On the other hand, the whole abelian I/(l) em x U(1)y gauge symmetry is restored 
in the core of a hybrid vortex (see Figs. 3 and 4). 

Recently, there has been a considerable interest in systems with coexisting order parameters (such as high T c 
superconductors) in condensed matter 0. Generating vector condensates is a very natural way of creating such 
systems (for example, in the present model, electric superconductivity coexists with spontaneous rotational symmetry 
breaking). This possibility deserves further study. 
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